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ABSTRACT 


5 “Pluctuations in the geomagnetic field were measured by 


- three orthogonally mounted coil sensors at two land sites 


separated by 40 km. Computer generated voltage vs time and 
RB magnetic field vs time plots failed to reveal the presence 
a of dominant micropulsations. A coherence study between the 
- two sites revealed coherence values of 0.6 - 0.8 in the 

- frequency range 0.04 - 0.6 Hz. This is compared to a 

~ 


coherence study completed at the Naval Air Development 
5 Center, Warminster, Pennsylvania, between land sites 24.8 km 


\ 


2 apart. The NADC coherence values are lower (0.3 - 0.6). 
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I. INTRODUCTION 


This thesis is part of an ongoing effort at the Naval 
Postgraduate School to analyze ULF geomagnetic noise and 
micropulsations. These variations in the geomagnetic field 
are of interest both from a geophysical and a military view- 
point. Applications of interest to the Navy are in the areas 
of magnetic detection of submarines, mine warfare and com- 
munications systems. 

The specific objectives of this study were to install i 
and operate a simultaneous data collection system at two 
separated land sites, to modify and adapt previously de- 


veloped software for data analysis and to obtain spectral 


coherences between the two sites for background noise and/or 
micropulsations. 

A coherence study of background noise with background 
noise, of micropulsation with micropulsation and of back- 
ground noise with micropulsation between the two sites 
should further the understanding of the types and extent 
of the sources that produce these fluctuations. 

The data collection sites were separated by a distance 
of 40 km (see Appendix A). One site was at La Mesa 
Village, near the Naval Postgraduate School campus, while 


the other was at the Chew's Ridge fire lookout. The latter 


was chosen for its remoteness from the local power grid. 


II. BACKGROUND 


A. MICROPULSATIONS 

The frequency spectrum of the geomagnetic field observed 
on or near the earth's surface has a number of well defined 
peaks, corresponding to categories of regular geomagnetic 
micropulsations, as shown in Figure 2.1. These micro- 
pulsations are designated as Pcl, Pc2, ... Pc5. 

Another category of micropulsations encountered is 
irregular pulsations. Unlike regular Pe micropulsations, 
which have relatively well defined frequencies, the Pi micro- 
pulsation consists of a spectral band of noise. 

The source of these micropulsations appears to be 
magneto-hydrodynamic resonances in the earth's magnetosphere 
(Pc2 ~ Pc5), ion cyclotron wave-particle interaction in the 
Magnetosphere (Pcl) and ionospheric currents perturbed by 
conductivity variations (Pi). References 1 and 2 give 
more detailed explanations of these mechanisms. Micro- 
pulsations are classified as follows: 

1. Pel: (0.2 - 5 Hz frequency) 

Known as "Pearls", these micropulsations are gener- 
ated by the cyclotron instability of energetic protons. 
They have been positively correlated with solar disturbances 


and occur during daylight hours in the auroral zone and 
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Figure 2.1 Field Strength of Micropulsations. 
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during night and early morning hours in the midlatitudes. 
Typical amplitudes are 0.05 - 0.1 nanotesla. 
2. Pe2: (0.1 ~ 0.2 Hz frequency) 

This is a diurnal phenomenon that shows some 
positive correlations with solar activity and the seasons. 
They usually decrease in their period as magnetic activity 
increases. Their average amplitude is 0.1 - 1 nanotesla. 

3. Pe3: (0.022 - 0.1 Hz frequency) 

These are similar to Pc2 pulsations except for the 
frequency range. 

4. Ped: (6.7 - 22 mHz frequency) 

Sunspot activity appears to have an effect on Pc4 
pulsations. Their frequency varies with the season and 
they have an average amplitude of 5 - 10 nanotesla. 

5. Pce5: (1.7 - 6.7 mHz frequency) 

These large scale pulsations occur during morning 
and evening with amplitudes of 10 - 100 nanotesla. Their 
duration shows a strong geomagnetic latitude dependence. 

6. Pil: (0.025 - 1 Hz frequency) 

These pulsations usually occur at night and early 
morning and vary in intensity from 0.01 - 0.1 nanotesla. 
They demonstrate a positive correlation with auroral 


disturbances. 


7. Pi2: (6.7 - 25 mHz frequency) 


The amplitude of these pulsations ranges from 1 - 5 


nanotesla. They usually occur during early morning hours 
but may continue throughout the night. The frequencies of 
these pulsations increase with increasing magnetic activity. 
Geomagnetic micropulsations can be distinguished from 
the general noise background of the geomagnetic field. The 
micropulsation events rise out of the ever present back- 
ground activity, reach an amplitude that can be large in 
comparison to the background level, and then finally disap- 
pear into the background. The Pc4 and PcS pulsations can 
last several hours. However, the Pcl - Pc3 and Pi micro- 
pulsations have a maximum duration of approximately one 


hour but may last only a few minutes. 


B. GEOMAGNETIC BACKGROUND NOISE 

It has been speculated that the primary source of the 
geomagnetic background noise is fluctuations in the inter- 
Planetary magnetic field [Ref. 3]. If so, a source of such 
large spatial extent implies that the amplitude of the 
background noise may be less variable over the surface of 
the earth than the more locally generated micropulsations, 


and one could expect considerable spatial coherency of the 


background noise over the earth's surface. 
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we David and Heirtzler [Ref. 4] studied the coherence of 
(. geomagnetic variations between two stations up to 550 km 
Soy apart. The geomagnetic variations were separated into a 


background noise component and a micropulsation component. 


. When two different micropulsation types occurred simul- 
“ taneously, they were found to be incoherent with one another 
> 
te and with the background noise. It would thus appear that 
2S 
= independent generation mechanisms exist for the background 
rs noise component and for micropulsations of different types. 
Co Also, the background component showed association with the 
he 
ad solar quiet day magnetic variation (Sq). In particular, 
ae the spectrum amplitude of the background component increased 
ke . 
oe as the strength of Sq increased. 
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III. DATA COLLECTION SYSTEM 


A. EQUIPMENT DESCRIPTION 
The system used at both the Chew's Ridge and La Mesa 
Village sites is shown in Figure 3.1. The major components 
are: | 
(1) Coil sensors 
(2) Preamplifiers 
(3) Signal conditioner 
(4) Pulse Code Modulation (PCM) encoder 
(5) WWV radio receiver 
(6) Tape recorder 
‘ (7) Power source 
For a geographical description of the two sites, see 
Appendix A. 
1. Coil Sensors 
Each coil is continuously wound with 5460 turns of 
18 gauge copper magnet wire. It has an internal resistance 
of 9.31 Henries. At each site, the three coils were 
mounted orthogonally, with the x coil oriented towards 
magnetic north, the y coil towards magnetic east and the 
z coil vertically downwards. 
2. Preamplifiers 
The preamplifiers are model 13-10A low noise 
amplifiers manufactured by Dr. Allen Phillips of SRI 
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a International. The overall power gain is 60 dB for inputs 
( less than 2.5 mv. A final stage low pass filter which 

Rn provides a sharp cutoff at 20 Hz is provided. Each pre- 
ie amplifier has a DC offset potentiometer which must be ad- 


ma justed to provide the correct zero-level at the output. 
3. Signal Conditioner 
The signal conditioner receives the analog signals 


from the preamplifiers, amplifies them by 30 dB and limits 


we them to an amplitude of 7.5 volts. 
re 4. PCM Encoder 
o The pulse code modulation (PCM) was designed and 
E manufactured by Dr. Robert Lowe of Lowecom Incorporated and 
= of the Scripps Institute of Oceanography in La Jolla, 
nea ° California. The encoder features 15 channel analog input 
we capability with selectable sample rates of an samples per 
oe second, where n is an integer value of 3 to 7. For the 
we purposes of this thesis only 3 of the input channels were 
iby utilized (1, 2 and 3 for the x, y and z channels, re- 
= spectively), and a sample rate of 64 samples per second 
a was chosen to ensure adequate measurement of frequencies 
a below 0.1 Hz. The encoder samples the analog signal from 
. channel at a rate of 64 Hz and assigns a pulse coded 
ne word with a decimal value between 0 and 4096 to 
es Guc.. sample, corresponding to an amplitude of -5 to +5 
mS volts. The output data is organized into frames, each 
mee 15 


OT 1385 


* 


READ, | 


“et 
ry 


PY 


& 
Lan 
. 


frame headed by a synch code word which is followed sequen- 


tially by the pulse coded samples from input channels 1 
through 15. The synch code word is a pulse coded digital 
word with a decimal value of between 0 and 4096. This word 
is preselected and hardwired on the encoder circuit board. 
Reference F explains the PCM system in more detail. 
5. WWV Radio Receiver 
In order to ensure that the data from the two sites 
was analyzed simultaneously, an R-1051 B/URR radio receiver 
was used to monitor the WWV Universal Time broadcast at 
20 MHz at each site. The broadcast gives the Universal Time 
at each minute by voice with each second marked by a tone. 
6. Tape Recorder 
Hewlett-Packard HP3964A/3968A tape recorders were 
used to record the PCM data and WWV broadcast on analog 
Magnetic tape. The output from the PCM encoder was recorded 
on a direct channel (100 - 16000 Hz frequency response) and 
the WWV time signal was recorded on an FM channel. 
7. Power Source 
At the Chew's Ridge site the power source used was a 
3500 watt, 60 Hz, 120 volt, gasoline powered, portable 
generator. The separation between the sensor coils and in- 
strumentation was about 100 feet, and between the sensor 
coils and portable generator approximately 250 feet. 
Commercial 60 Hz power was available at the La Mesa Village 


site. 
16 


The preamplifiers, signal conditioners and PCM encoders 
were powered by rechargeable 18 amp-hour batteries (plus 


and minus 12 volts and ground). 


B. PCM TO DIGITAL CONVERSION 

The PCM data recorded on analog tape is played back into 
a PCM decoder which converts it to digital data. The digital 
data is recorded in 9-track, 800 bits per inch computer 
tape for subsequent analysis on the IBM 3033 mainframe 
computer. Appendix B contains a step-by-step procedure 
for the decoding process. 

By listening to the FM channel carrying the WWV time 


signal over a speaker, the point on the analog PCM tape 


where it is desired to begin and stop the decoding process 
may be precisely determined. In this manner it is pos- 
sible to obtain time synchronized digital computer tapes of 


data from the two sites. 


IV. COMPUTER SOFTWARE 


The computer programs used to analyze the data are 
written in Fortran IV programming language and are briefly 
discussed below. These programs are listed in Appendices 


Cc - G. 


A. PROGRAM VOLTR 

The VOLTR program reads data from a digital computer 
tape and generates a voltage vs time plot for each orthogonal 
axis. The data is read from the tape in blocks of 8192 
frames (128 seconds) by the subroutine RD. This data, 
which is in integer form between 0 to 4096, is then normal- 
ized to represent voltages between + 5 volts. The amount of 
data plotted is an integer increment of 128 seconds, the 


integer being from 1 to 8 and specified by the user. 


B. PROGRAM VODIG 

This program applies a 144 point double running average 
and a .04 - 0.6 Hz digital filter to the rough voltage and 
generates filtered voltage vs time plots for each axis. 
The digital filter simulates the pass band of an AN/SQ-81 
Magnetometer and was developed by Mike Huete of the Naval 
Postgraduate School. Reference 5 explains the filter in 
detail. The double running average smooths out any large 


noise "spikes" that may cause an unnatural oscillatory 


18 


response in the digital filters. It also acts as a low pass 


nl filter, removing frequencies greater than approximately 1 Hz. 


C. MASS STORAGE PROGRAM 

In order to compare simultaneous data from two different 
computer tapes, the data is read from one tape (La Mesa 
Village), normalized to voltage values and stored in the 
IBM 3033 Mass Storage System, where it is available for 


future recall. 


D. PROGRAM MAGFLD 


Tee! 


This program generates magnetic field vs time plots. 
The digital data is read from the computer tape and normal- 
ized to voltage values. A Fourier transform is performed 
: on the data to enter frequency space. At this point the 

system transfer function, which converts the data from 

aay voltage to magnetic field values, is applied. References 
6 and g detail the procedures used to determine the trans- 
fer function for each coil sensor-amplifier subsystem. 
After the transfer function has been applied, a second 
Fourier transform is performed to return the data to time 
space. A 144 point double running average is then applied 


- to the magnetic field data to remove frequencies above about 


4 A 1 Hz. 
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B. PROGRAM COHER 


This program calculates the spectral coherence of the 
total field between the two sites and the power spectral 
densities of the total field at each site. The La Mesa 
Village data previously stored in the Mass Storage System 
is recalled, the corresponding Chew's Ridge data is read 
from a computer tape, and the two data sets manipulated 
simultaneously. 

Referring to Figure 4.1, the total field was calculated 
as 


total field = xcosé, + zsinO, 


where 94 is the magnetic dip angle, which in the Monterey 


area is 60°. 


> TF 


Figure 4.1 Total Field Diagram 


The coher_nce between two signals a(t) and b(t) is 


a(t) eb(t) aati A(F) B*(F) _ 
YValt)ealt) wb(t)eob(t) Va(F)A*(F) YB(F) B* (F) 


where 'o)' indicates the correlation separation, '*' indicates 


coherence = 


20 
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the complex conjugate is taken and A(F) and B(F) are the 
Fourier transforms of a(t) and b(t) , respectively. 

In the program, a Fourier transform of the data into 
frequency space was performed and an average of 20 blocks 
of data (128 seconds per block) was taken to obtain the 


final coherence values. 
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V. TESTING OF PCM SYSTEM AND SOFTWARE 


In order to ensure that the PCM svstem and VOLTR program 
faithfully reproduced the input signals, sinusoidai, tri-~ 
angular and square waves were input to the PCM encoder by a 
Wavetek signal generator, as shown in Figure 5.1, and the 
PCM signal recorded on analog tape. The signal generator 
output was also monitored by a chart recorder and voltmeter. 
The analog tape was then decoded and voltage vs time plots 
were generated by the VOLTR program. 

Table 5.1 shows the relationship between the amplitude 
of the signal generator output and the amplitude of the 
VOLTR plots (Figure 5.2) at various frequencies for the 
sinusoidal signal. 

As can be seen from Table 5.1, the error between the 
chart record and the computer plot is less than two percent. 
Similar results were obtained for the triangular and square 
wave on all three channels. 

Extensive testing of the Mass Storage program, and the 
digital filter algorithm employed in the VODIG program, is 
documented in References 6 and 7, respectively. 

The program COHER was tested by analyzing a section 
of data against itself. Data from a computer tape was 
read into the mass storage system by the Mass Storage pro- 


gram. The same section of data was read from the tape by 
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a 

ee the COHER program and analyzed with the data recalled from 
[ Mass storage. The COHER program generated a coherence vs 
= | frequency of 1 as expected (Figure 5.3). 

= Reference 8 mentioned the presence of "cross-talk" 

me between the channels of the PCM encoder. This was noticed 
es on computer generated plots on a channel whose input jack 
BS was left open while making 2asurements at a field site. 
— To test for "cross-talk", a signal from the Wavetek signal 
xs generator was fed into all three channels of the PCM en- 
= coder, as in Figure 5.1. One of the channels was dis- 

if connected from the Wavetek and the input jack left open. 
= Then the input jack was grounded, and then finally the 

ce Wavetek signal was reconnected. Figure 5.4 is a rough 

= voltage vs time plot for the sequence. It can be seen that 
SS while the input jack was open a signa} did appear on the 
= channel but disappeared while the input jack was grounded. 
an The "cross-talk" mentioned in Reference 8 was actually the 
ie open input jack acting as a "pickup" antenna. 
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VI. EXPERIMENTAL RESULTS 


Data was taken on 4 August 1983 between 1300 and 1845 
local time. Beginning and ending the recording of the 
analog tapes was coordinated between the two sites over 
PRC-77 radios. Since it proved difficult to communicate 
between the two sites directly (because of intervening 
hills), a person using a radio with a large whip antenna 
at the Naval Postgraduate School directed the simultaneous 
starting and ending of data recording at the two sites. 

On the voltage and magnetic field plots, the units 
labeled on the vertical scale are arbitrary and only the 


peak-to-peak variations should be considered. 


A. ROUGH VOLTAGE PLOTS 

Figures 6.1 - 6.6 show the rough voltage plots for 
the La Mesa Village site. These signals are totally 
obscured by 60 Hz noise. Figures 6.7 - 6.12 show the 
Chew's Ridge rough voltage plots. Here the 60 Hz noise is 
a site as remote as Chew's Ridge (to escape the 60 Hz 


power grid) is thus justified. 


B. FILTERED VOLTAGE PLOTS 
Figures 6.13 - 6.30 show typical filtered voltage vs 


time plots for both sites. Visual inspection failed to 
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reveal the presence of any large amplitude micropulsations 


or of any clear one-to-one correspondence in simultaneous 


sections of data. 


C. MAGNETIC FIELD PLOTS 
Figures 6.31 - 6.45 show typical magnetic field vs 
time plots for both sites. Magnetic field variations at 
the La Mesa site are approximately one nanotesla; variations 
at the Chew's Ridge site are slightly greater, 2 - 4 


nanotesla. 


D. COHERENCE PLOTS 

Figures 6.46 - 6.57 show coherence vs frequency plots 
for individual axes and for the total field. The coherence 
generally has values between 0.6 - 0.8 indicating a moder- 
ate degree of commonality in the geomagnetic variations at 
the two sites. 

These coherence plots can be compared with coherence 
vs frequency plots generated from background geomagnetic 
variation data taken at the Naval Air Development Center 
in 1979. The separation between the NADC data collection 
sites was 24.8 km. Figures 6.58 - 6.60 show these plots. 
In general, the coherence values from the NADC data are 
less than the coherence values found in our measurements. 


However, the amplitudes of geomagnetic variations are 
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probably influenced by factors such as the state of the 
ionosphere and magnetosphere and the stage of the solar 


cycle. 


The NADC data was averaged over a period of two hours 


while our data was averaged over a period of 40 minutes. 
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VII. CONCLUSIONS AND RECOMMENDATIONS 


A coherence of 0.6 - 0.8 in the background component 
of the geomagnetic field between 0.04 and 0.6 Hz was 
established. A coherence in this range should be considered 
only moderate. A lack of high coherence (0.9 - 1) indicates 
that the variations in the background field observed at the 
earth's surface at the two sites are not produced directly 
by the same source. However, the variations are clearly 
not random in nature. The moderate coherences found suggest 
that the source mechanisms for the background component in 
the geomagnetic field are complex and involve mechanisms in 
addition to or intermediate to simple fluctuations in the 
interplanetary magnetic field. 

A discernable micropulsation was not recorded during 
the five hours of data taken. It is recommended that 
additional data be taken at the two sites in the hope of 
performing a coherence study on the micropulsation com- 
ponent of the geomagnetic field. It is also recommended that 
data be taken at additional sites of greater separation 
(100 km or more) in order to investigate the degree of 


coherence with distance. 
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APPENDIX A 


SITE DESCRIPTION 


The Chew's Ridge fire lookout is located 40 km south- 
east of the Naval Postgraduate School and at an altitude 
of approximately 3900 feet above sea level. It was chosen 
for its remoteness from the local power grid. Since the site 
is within the Los Padres National Forest, permission to 
collect data there had to be obtained from the National 
Forest Service. A dirt road provides easy access to the 
site for the transporting of equipment. The Monterey 
Institute for Research in Astronomy is currently constructing 
an observatory approximately one half mile from the lookout. 
What affect its presence will have on the suitability of 
the fire lookout for future data collection is not currently 
known. 

Initial attempts to transmit the PCM data via a 170 MHz 
carrier wave from this site to the school proved impossible 
due to the relatively low transmission power used (3 watts) , 


less than ideal line of sight, less than ideal antenna. 


93 


bd hel bet het ll nn 
. _ . nea. ‘ 


welt, 
M@ * *.f. 


Figure A.l Geographical Area of Data Collection. 
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APPENDIX B 


PCM DECODING PROCEDURE 


Several electronic components are utilized to decode the 
PCM data to digital form which is ultimately stored on nine 
track, 800 bits per inch digital tape. This data pro- 
cessing system is shown in Figure B.1l. Central control of 
this process is accomplished with a Hewlett Packard 9845A 
computer utilizing an operator interactive program, "PCMPROG". 
After execution of the program, the computer requests entry 
of specific function control parameters into the computer 
and other equipment. These inputs are used to control synch- 


ronization of equipment start, digital tape drive speed, 


decode rate, decode time and synchronization code word 


entry into the decoder. The PCM encoded data is fed into 
the system from the HP 3964A/3968A tape recorder previously 
used to record the data. The decoding of the PCM data is 
accomplished with a Marine Profiles, Incorporated Model 

319 PCM decoder. A Monsanto AM-6419/USM-368 oscilloscopes 
are used to display the PCM data. A Kennedy Model 9800 
digital recorder and computer interface are employed to 


store the digital data on the nine track digital tape. 
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Figure B.1 Decoding System Data Flow and Control. 
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1. Decoding Procedure 


a. Connect a coaxial cable between the output 
channel desired for decoding of the HP 3964A/3984A tape 
recorder and channel three ofthe Model 319 PCM decoder. 

b. Energize the HP 3964A/3984A tape recorder, the 
Monsanto AM-6419/USM-368 oscilloscope, the Kennedy inter- 
face and the AANDERAA tape transfer interface, the Model 319 
PCM decoder and the HP 9845A computer. Also, on the PCM 
Model 319 decoder place the fan toggle switch in the UP 
position and the AC/DC/OFF switch in the AC position. 

c. Place into the right hand side tape reader of 
the HP 9845A computer the program named "PCM PROG". Type 
the command "MASS STORAGE IS":T15" and press EXECUTE. 

Then type the command GET "MT" and press EXECUTE. 
dad. On the PCM decoder place the following functions 
to the listed positions: 
SOURCE - 3 
SAMPLE RATE - 64 (for 3 3/4 recorder speed) 
INVERT/NORMAL - NORMAL 
OUTPUT/SAMPLE RATE - 0 
RECORDS/FILE - INFINITY 
SYNC CODE - 000 
e. Press RUN on the computer, ignore the computer's 


response "enter Y to skip tape init" and press CONTINUE. 
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E: tid it 
< 

{ 

= | f. The computer now indicates "load tape" into 

a the Kennedy unit and "put on line". ees this energize 
= the Kennedy unit, place a write ring on the digital tape 
s and load the tape according to the diagram located on 

= the inside of the unit's door, press the LOAD button and 
s the ON-LINE button located on the front of the unit. 

a g. The computer now indicates "enter synch code". 
. Type into the computer 3658 (Chew's Ridge tapes) or 3155 
= (La Mesa Village tapes) and press CONTINUE. 

3 h. Enter transfer time in minutes and seconds into 
ia the computer. For example 30 minutes and 50 seconds would 
= be typed in as 30,50. Most analog tapes lasted 45 - 50 

- minutes. After this is done, press CONTINUE. 

a : i. Push the STOP switch on the PCM decoder and 

% press CONTINUE on the computer. 

yi j. Push the PLAY button on the Hewlett Packard 

; tape recorder, listening to the WWV time signal over a 

~ speaker or headphones. Push the START switch on the PCM 
- decoder to begin the decoding process at a chosen time, 

“a using the second "ticks" of the time signal as a count- 

vs down. The decoding of the corresponding analog tape from 
S the other recording site must be started at precisely the 
rq ; same time. 

< 
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k. To end the data transfer early, push the KO 
button on the computer. If this option is selected, "T" 
must be entered on the computer to write end of file on 


the digital tape. 


1. “End of run” will be indicated on the computer 


CRT. Deenergize the equipment. 
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